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ABSTRACT 



The finite-temperature renormalization group is formulated via the Wilson-Kadanoff 
blocking transformation. Momentum modes and the Matsubara frequencies are coupled by 
constraints from a smearing function which plays the role of an infrared cutoff regulator. 
Using the scalar Xcp 4 theory as an example, we consider four general types of smearing 
functions and show that, to zeroth-order in the derivative expansion, they yield qualita- 
tively the same temperature dependence of the running constants and the same critical 
exponents within numerical accuracy. 



CCU-TH-98-01 March 1998 



' electronic address: senben@phy.ccu.edu. tw 
* electronic address: strickland.41@osu.edu 



1 



I. INTRODUCTION 



Finite-temperature renormalization group (RG) techniques provide a powerful means 
of probing the physical behavior of a quantum system interacting with an external heat 
bath. The presence of the additional length scale (3 = T -1 , the inverse of the temperature, 
makes the theory more complicated since there are now two types of fluctuations, quantum 
and thermal, which compete with one another. In order to apply the RG framework to 
explore issues such as the QCD deconfinement phase transition and the evolution of the 
early universe, one must first understand the subtle interplay between these fluctuations. 

To take into account both thermal and quantum fluctuations, various finite- 
temperature RG prescriptions have been proposed [1] , [2]. These RG schemes also 
provide an improvement beyond perturbative approximations, which are known to be 
plagued by infrared (IR) singularities in the high-temperature limit as well as in the vicin- 
ity of criticality. In the present work, we revisit the studies of Ref. [1] which were based 
on the Wilson-Kadanoff blocking transformation [3] performed on a manifold S 1 x R d in 
the imaginary-time formalism. By introducing a regulating smearing function p^(x, t x ), 
which simulates the coarse graining in momentum space, the effective degrees of freedom 
of the theory are now characterized by the blocked fields, 4>n,k, 

(pn,k(^ T ^) = j o dT V f P^l(^-y^x-T y )(j)(y,T y ) , (1.1) 

where <f>(x) represents the original field variables. While n is a discrete cutoff for the 
Matsubara frequencies u n = 2nn/(3, the scale k may be regarded as a generalized IR 
cutoff for a combination of both |p| and u n . Varying k infinitesimally then allows us 
to generate a nonlinear differential RG equation for the blocked action Sp^\<&\, which is 
defined by 

e -s fi , k [*M]= f D[(j>} n<S(</> fc (x)-$(x))e- s ^. (1-2) 

J periodic x 

The blocked action provides a smooth interpolation between the bare action Sa[$] and 
the quantum effective action which generates the one-particle-irreducible graphs as k is 
lowered from the ultraviolet (UV) cutoff, A, to zero. We have taken the field average of a 
given block <fr(x) to coincide with the slowly varying background. Since the flow equation 
for Sp t k[&] is too complicated to be solved exactly, we make an expansion in powers of 
derivatives: 

Sf,,k[*] = f Q *r / x {^W) 2 + ^|M(V$) 2 + U,,m + 0(d 4 ) j . (1.3) 

The theory is then parameterized by a set of coupled nonlinear partial differential RG 
equations for the wavefunction renormalization constant, Zp^{f&), the blocked potential, 
Up,k(®) and higher-order contributions. The behavior of the theory at any arbitrary k and 
T can then be studied by analytically or numerically solving these coupled equations. We 
shall designate this type of RG prescription as the momentum RG, or MRG. 
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Notice that the functional form of p^ £(x, r x ) in MRG is not unique; different choices 
will change the way in which the irrelevant, short-distance degrees of freedom are eliminated 
in course of blocking. In the finite-temperature setting, how the two types of modes, u n 
in S 1 and |p| in R d , are block averaged is controlled by p^(x, r x ). In Fourier space, 

p~ k(p,u; n ) modifies the bare propagator as 

1 PnjiP^n) 

Thus, one may visualize the following possibilities: (i) Blocking is performed in only one of 
the two subspaces, i.e., with Pn(r x ), or (x); (ii) both u n and |p| are blocked, but with 

two smearing functions Pn(r x ) and pj^(x) that are independent of one another; and (iii) a 
constraint exists in the smearing function for u n and |p| which allows them to be coupled 
and eliminated in some coherent manner. The trivial case of no blocking simply implies 
the absence of any RG transformation. Scenario (iii) is the most interesting case because 
it allows us to examine the interplay between the quantum and thermal fluctuations in the 
two spaces. In this paper smearing functions corresponding to the above categories are 
constructed explicitly. In particular, we study, in Fourier space, p\f\p) = Q(k — |p|) for 
category (i), p^ k (p,oj n ) = 0(k-\p\)Q(h-\n\) for (ii), and p k d \p,oj n ) = Q(k-^/u% + |p| 2 ) 

and p k (p, u n ) = Q(k — \ui n \ — |p|) for (iii). These smearing functions differ only in their 
treatment of the coupling between |p| and the n ^ sector. Thus, by comparing the 
numerical solutions of the corresponding RG equations, one can gain insight into the 
importance of the feedback between u n and |p|. 

To derive the MRG equations, we lower k infinitesimally from k to k — Ak, where Ak is 
the characteristic thickness of the RG shell that contains the modes to be eliminated. The 
ratio k = Ak/k represents the fraction of the modes that are eliminated in each blocking 
step. When the smearing function involves n explicitly, the resulting RG becomes piece- 
wise continuous due to the discreteness in the Matsubara frequencies. However, we shall 
find that in all the cases considered, qualitatively the same results are obtained for the 
temperature dependence of the running parameters. The quantitative difference stems 
from the boundary conditions on the RG flow when the Matsubara modes are blocked. 

We also show that all four methods yield the same critical exponents to within 0.5% 
when Uf3 y k(&) is expanded in a power series in $ near criticality. Thus, it may be concluded 
that the feedbacks between |p| and the n^O sector can be treated perturbatively in the 
high T regime. Although scaling patterns in the IR regime differ in all four cases, direct 
correspondence may nevertheless be established and equivalence may be recovered in some 
cases. 

We emphasize that k must always remain small in order to trace out the true RG 
evolution. A large k causes poor tracking of the RG trajectory. To illustrate this point, 
we examine the high-temperature behavior of the quartic coupling constant of the scalar 
A</> 4 theory in the symmetric phase. It has been shown that A^ = A / 3 i / £ =o decreases with 
increasing T but approaches a constant using p k d \p) = Q(k — ]p|) [4]. However, if one 
employs the "Matsubara shell assumption" with Ak = 2n//3 , i.e., the thickness of the RG 
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(1.4) 



shell is chosen to be equal to the spacing between two adjacent Matsubara frequencies, an 
incorrect high-temperature behavior will be obtained. The result is due to the fact that 
as Ak becomes increasing large in the high T limit, k is no longer small and the flow then 
deviates significantly from its true trajectory. 

The organization of the paper is as follows: In Sec. II we give four examples of the 
MRG smearing functions and describe the details of how the corresponding RG evolution 
equations are generated. We compare and contrast the manner in which u n and |p| are 
coupled. Numerical solutions for the MRG equations are presented in Sec. III. We demon- 
strate once more how decreases with increasing T and may be approximated by 
a constant in the infinite-temperature limit. Results extracted from the Matsubara shell 
assumption are also included to illustrate its failure in the high-temperature limit. Sec. IV 
is reserved for summary and discussions. In Appendix A we present a complementary RG 
scheme, the temperature RG (TRG), which is formulated using T as the running param- 
eter. Its consistency with MRG is demonstrated. In Appendix B the connection between 
the the finite-temperature behavior of the running parameters and the fixed points is es- 
tablished. In particular, we identify all the fixed points accessible by both the MRG and 
the TRG schemes and show that the existence of an infinite-temperature Gaussian fixed 
point rules out the possibility of \p,k=o increasing with T. 

II. MRG SMEARING FUNCTIONS 

As described in the Introduction, under blocking transformations an IR cutoff k is 
introduced into the theory by defining the coarse-grained fields (1.1) via a smearing function 
p~/ k {x., t x ). Here h differentiates between the high and the low Matsubara frequencies, 
and k provides a separation between the slowly varying background fields and the fast- 
fluctuating modes which are to be integrated over. Below we consider four smearing 
functions, each of which has a different coupling between u n and |p|. 

(1) pi d) (p)=e(fc-|p|): 

The simplest blocking is to choose a sharp cutoff in the R d submanifold, leaving S 1 
unconstrained. In position space, this smearing function has the form 

Pf(x) = f e'P~ (2.1) 
J\ P \<k 

which approaches 5 d (x) as k — > oo. This limit corresponds to a completely "unblocked" 
bare system characterized by the original unrenormalized fields. The blocked fields take 
on the form [4] 

^(x,r x )= /pi d) (x-y)0(y,r x ) = i £ / e"^—- ^(p). (2.2) 

We see that the effect of this smearing is to provide an integration over |p| > k for all n. 
However, with this choice of p k (x) , no averaging is done in the r x direction. 
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Let us first consider the simple, independent-mode approximation where all the modes 
greater than k are eliminated at once with no feedback being taken into account. The 
perturbative one-loop correction to the blocked potential reads: 



2/3 



E 



In 



+ P 2 4 



^ J |/Vp 2 + V»(*) + 2 ln[l - e -^P a +v»(*)] | + 



(2.3) 



where the prime notation in the integration implies k < |p| < A. The wavefunction 
renormalization constant has been set to unity for simplicity. A direct differentiation of 
the above expression with respect to k then gives 

k du^m = _s^_L^ k2+Vfl ^ )+2 ln [i _ e -ev^W)} \ , (2 .4) 

dk 2p L J 



where = 2/(47r) d / 2 r(<i/2). Unfortunately the equation has a serious drawback in that it 
only accounts for the one-loop contribution and neglects the important feedback from the 
higher to the lower modes whose energy scales, in principle, may differ by several orders 
of magnitude. Since the effective degrees of freedom of the theory are not well tracked, 
it will fail to capture the relevant physics in the high temperature limit as well as in the 
neighborhood of the phase transitions. 

The aim of MRG is precisely to take into consideration the important higher loop 
corrections and the couplings between the modes. In Figure 1 we depict how RG is carried 
out as an improvement beyond the independent- mode approximation. At k = A, the 
system is unblocked and no modes are eliminated. As the system undergoes blocking with 
the IR cutoff being lowered from k to k — Afc, the modes contained inside the RG shell of 
thickness Ak are integrated out. The elimination of these fast-fluctuating modes results in 
a more complicated blocked potential Up^-Aki®) which contains more operators than the 
previous one, Up^{Q). The evolution is complete when k finally reaches zero. By lowering 
k infinitesimally, the RG equation reads: 



k 



dUp, k {$) __ S d k d 



dk 



2/3 ^ 



E "» 



2 + k 2 + U" 



p,k 



(2.5) 



which, upon summing over n, yields 



k 



dUp, k (®) 
dk 



S d k d 



(2.6) 



The functional forms of Eqs. (2.4) and (2.6) are remarkably similar, except for the sub- 
stitution on the right-hand-side of the latter by the scale-dependent counterpart, i.e., 
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V"(3>) — > Up k (<&). Though this "dressing" seems straight- forward, it has been demon- 
strated to have far-reaching consequences, particularly in the study of phase transitions 
and critical phenomena. By solving the equations numerically subject to the initial con- 
dition £//3 =00) a(<I>) = Va(<E>), where Va($) is the bare potential defined at the UV cutoff 
A, the severe IR divergences encountered in the finite-loop calculations are completely re- 
moved. In addition, critical exponents can be measured to a remarkable precision through 
numerical solutions of the exact renormalization group flow equations [5]. 



RG Shell □ Eliminated modes 




Ak 



Figure 1. Schematic diagram of blocking Method 1. 



It is useful to compare Eq. (2.6) with the zero-temperature result for D dimensions, 
obtained using a sharp cutoff for .D-dimensional momentum integration [6]: 



k 



dk 



Snk D 



In 



k 2 + U£W 



(2.7) 



In the regime where (3 
and one has 



0, the contribution from the first term in (2.6) can be neglected 



k 



dUf3, k ($) _ S d k d 
dk 2/3 



In 



k 2 + U'^ k m 



(2.8) 



which is simply the contribution from the n = mode, as can be seen from (2.5). Thus, 
we see clearly that the n = sector alone is sufficient to account for the <i-dimensional 
characteristics. Upon rescaling of $ — > /3 _1 / 2 $ and Up t k ~^ /3~ 1 Uk : Eq. (2.8) coincides 
with Eq. (2.7) for D = d. 

On the other hand, in the low T regime where (3 — > oo (holding k fixed), a com- 
plete summation over n is necessary for recovering the proper zero-temperature (d + 1)- 
dimensional feature of the theory. The expected RG equation is that of Eq. (2.7) with 
D = d + 1. However, in this limit Eq. (2.6) gives 



k 



dU k ($) S d k c 
dk 2 



(2.9) 



6 



instead. A careful comparison of the RG flow governed by Eqs. (2.7) and (2.9) shows that 
they lead to the same renormalized values for the vertex functions at k = 0, but differ in 
the scalings patterns for k <C where pn is the renormalized mass of the theory [4]. 
How can we improve our RG to ensure that the desirable IR scaling limits below the mass 
gap are attained? To do so, one must carefully take into account the couplings between 
uj u and |p|. From Figure 1, it is clear that Method 1 results in a complete summation over 
all Matsubara modes within a given RG shell and ignores the couplings between modes 
having the same |p| but with different n's. How do we take into consideration this coupling, 
and how important is this to physical applications? To answer these questions, we must 
examine other smearing functions. 

(2) p^ k (p,u n ) = e(k-\p\)e(n-\n\): 

As an improvement to account for the systematic feedback between certain d- 
dimensional momentum modes |p| and the Matsubara frequencies u n , we choose the smear- 
ing function to be p^\(p,u n ) = Q(k — |p|) Q(h — \n\), or 

pS(x,t,) = Pn(T x )p { k d) (*) = I E / e- iKlT *- p ' x) , (2-10) 

which implies: 

fc (x,r x ) = iE / e-^-^Mp) • (2-11) 

P\n\<n J M< k 

In other words, this smearing function provides a summation over n > \n\ in S 1 and an 
integration over all |p| > k in R d . From 

1 oo 

n= — oo 

one can readily verify again that the limits k — > A and h — > oo correspond to p-^x, r x ) — > 
S d (x.)S(r x ) and the theory coincides with the original "unblocked" system. With such a 
step-by-step elimination, the couplings between different Matsubara frequencies are also 
taken into account [7]. 

It turns out that the bounds on n can actually be chosen in a more physically motivated 
manner. To do so, we notice that just as the high modes with |p| > k in R d are being 
integrated over, we should also sum over in S 1 all the high Matsubara frequencies with 
\oj n \ > k, or \n\ > [f3/2ir], where f3 = f3k and [j] is the greatest integer less or equal to j. A 
similar upper bound, (3A/2tt on cu n , or equivalently, n max = [@A/2n] may also be imposed. 
The blocking procedure of this scheme is depicted in Figure 2. 
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RG Shell 



Eliminated modes 




Figure 2. Schematic diagram of blocking Method 2. The thickness of the RG shell in the n-axis is 
zero for Afc^O (exaggerated here), but a Matsubara mode is eliminated after M continuous 
iterations where [N (3 Ak / 2tt\ = 1 . 

From the figure, we see that the RG shell on the uj n — |p| plane can be subdivided into 
three parts: one vertical strip having a thickness Ak containing modes with \n\ < [/3/2ir], 
and two horizontal strips with |p| < k. In the limit Ak — > 0, the contribution from the two 
horizontal strips vanishes, as long as Ak does not span a Matsubara mode. It is crucial 
not to make the Matsubara shell assumption with Ak = 2n / [3 such that all three volumes 
contribute. Such an assumption will lead to serious difficulties, particularly in the high 
temperature regime where Ak becomes so large that the treatment would be reduced to 
that of the independent-mode approximation. Thus, one always makes Ak small to ensure 
that the RG shell contains as few modes as possible. The subtraction of a Matsubara 
mode takes place only after every M continuous iterations of RG in the |p| space when 
MAk = 2tv/(3 and the condition [Af(3Ak/2ir] = 1 is reached. Similar blocking iterations 
are then repeated until we arrive at n = 0. 

The above prescription leads to 



Uf3,k-Ak - U/3 yk 



1 

w 

Sd 
2(3 



or, in the limit Ak — > 0, 



E / e(|p|-fc + Afc)e(fc-| P |)inLj+^ + c^ >fc (*) 



" <[&] 



s Jk-Ak J 



(2.13) 



k 



dk 



2(3 



[0/2n] 
= -[/3/2tt] 



(2.14) 
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We see that contrary to Eq. (2.6) in Method 1, the summation over n is now restricted 
to that contained in the shell, namely \n\ < [[3/2n]. This RG coincides with Eq. (2.8) as 
j3 — > 0. On the other hand, when j3 — > oo, one recovers (2.5) or (2.6) which contains both 
zero and finite-temperature sectors for d + 1 dimensions. However, the latter contributes 
only perturbatively [5]. 

The situation becomes slightly more complicated when one arrives at the step where 
a particular mode along the n-axis, say ±uj m (n = ±m), is contained in Ak and is to be 
eliminated. In this case, the RG equation becomes 

Uf,, k -* k {*)-U fl , k {*) = ± £ fle(\p\-k + Ak)Q(\n\<[^]) + 6 nt±m \ 

n= — oo P ^ ' 

x0(fc-|p|) \n[u> 2 n +p 2 + Ul k (<S>) 
(Ak) lnk+p 2 + ^ fc ($) 



Sjk d ~ 1 



2/3 



(2.15) 



l»l<[£]-i 



+ 



3d f , d-1 

Tk dpp 



In 



where the last term represents the contribution from ±u m . For d = 3 we can evaluate the 
last term analytically 



27T 2 Jo 

1 



6tt 2 



dp p 2 ln 
fc 3 ln 



UJ, r 



+ p 2 + ^ fc ($) 



k 2 + uj 2 m + U% 



+ 2[^ + ^] 3/2 tan" 1 (- 



k 



(2.16) 



+ 2k[u 2 m + Ul k ^)} 



The absence of the Ak factor in the expression prevents us from generating a differential 
equation for Up^ k (^). Nevertheless, Eq. (2.15) can be solved numerically. 

(3) Pk d \p^ n ) = Q(k - ^l + WY 

The above smearing function represents another class of finite-temperature blocking, 
and has been considered in Refs. [8] and [9]. In position space, it corresponds to 



p[ d \*,T x ) = - / e^—^) Q(k - y/u* + |p|2) . (2.17) 

This particular choice of p k d ^ (p, ui n ) leads to a mixing between the components in the 
subspaces 5" 1 and R d , and the Matsubara sum over n is now constrained by the value of 
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k, similar to Method 2 but in a more complicated manner. Figure 3 depicts how modes 
contained in the RG shell are eliminated as k is lowered. 

| | RG Shell O Eliminated modes 




Figure 3. Schematic diagram of blocking Method 3. 



If the RG shell does not contain a Matsubara mode along the n-axis, a nonlinear 
differential flow equation can be obtained as 



k 



dk 2/3 



where 



E / Ink + P 2 + Ug tk ($j\6(y/u,* + |p| 2 - k) 



S d k 2 
" 2/3 

S d k 2 
" 2/3 

S d k d 
" 2/3 



E / * 



d-i 



n/A 7 ' - ,.- 



2 _ ,.,2 
n 



In 



d-2 
2 



ln[^ 2 + ^($)]E( fc2 -^) 



</ d )(£) ln[fc 2 + ^ fc ($) 



9 M(« = l + 2^[l-(^) 2 



n=l 



d-2 



(2.18) 



(2.19) 



is a piece-wise continuous function of /3, as depicted in Figure 4 for d = 3. 

On the other hand, as in Method 2, an additional contribution representing the \m\- 
th Matsubara mode along the n-axis appears during their elimination, and the equation 



10 



reads: 



U/3,k-Ak($) -UfrkiQ) = 



oo 



2(3 * 



/ e(fe- ^ 

J O 



+ |P| 



) 



+ |P| 



k + Ak) 



+ 6 n ,± m e(s/w* + \p\*-u> m ) In u/J+p 2 + 



S,k d ~ 



2(3 



(Ak) gWtf) \n\k 2 + Ul k ($) 



Sd 



dp p 



d-i 



In 



(2.20) 

where g^ d \(3) is the same as g^ d \(3) in (2.19) except the n = ±[(3/2tt] contribution is 
excluded in the summation. Thus, we apply Eq. (2.20) whenever a particular \u m \ is 
contained in the RG shell. It could happen that the modes lie exactly at the boundary of 
the shell, i.e., k — Ak = u) m . In this case one may easily verify that Eq. (2.20) reduces to 
Eq. (2.18) and a differential flow equation is recovered. Although one could in principle 
fine tune Ak with respect to T in a way such that the points along the n-axis always lie 
on the boundary of the shells, we shall treat Ak and T as independent input parameters. 



15.0 



10.0 



to. 



0.0 



0.0 



50.0 
P 



100.0 



Figure 4. Plot of g (3) (0) for Method 3. 



The novel feature of employing this type of blocking is that for both (3 3> 1 and 
j3 1, one will obtain the expected zero-temperature result as given in Eq. (2.7), namely, 
D = d+ lfor(3^>l and D = d for (3 1. This can be readily seen by noting that in 
the high-temperature limit (3 — > (decreasing (3 while holding k fixed), g^ d \(3) — > 1 and 
Eq. (2.18) is reduced to Eq. (2.8). On the other hand, in the regime where (3 3> 1, the 
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summation in Eq. (2.19) can be replaced by an integral: 



g {d \(3) « 1 + 2 



i: 



dn 

d 

(4) 



,27m, 2 



d-2 
2 



d 3 4tt 2 n S, 



(2.21) 



2' 2' /3 2 



d+l 



(3 + 



In this limit, S 1 x R d is effectively R d+1 and the RG equation is reduced to that of (2.7) 
with D = d+l. 

(4) p{ d) (p,u n ) = e(k-\u n \-\p\): 

The above smearing function represents yet another method of blocking at finite tem- 
perature. The manner in which the modes are eliminated is illustrated in Figure 5. 

□ RG Shell □ Eliminated modes 




Figure 5. Schematic diagram of blocking Method 4. 



In analogy to Method 3, the components |p| and u n are also mixed due to the con- 
straint from the blocking function. By comparing the results obtained from these two 
schemes, one can gain insight into the effects of this mixing. The RG equations can be 
written as 



k 



dUp, k (<f>) _ k 



dk 



In 



u/J + p 2 + l/£ fc (*) < 9(n)8(uj n + |p| - k) 



+ <y„, <y(|p| -k) + 0(-n)6(-u n + |p| - k) 



K n=l 



27TO\ d - 1 1 

In 



> 



+ (k- u n f + Ul h <*)\ + ln[k 2 + U$ tk (*)] | , 

(2.22) 
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and 



k\ 



(3,k-Ak 



U, 



0,k 



Ak 



S d k d 
2/3 



[^/2tt]-1 , 
n=l 







wJ + (fc-u/„) 2 + Ctf fc (*) 



k—U> r 



+ ln[k 2 + Ul k {*)] | + J o dp p*- 1 ln[u& +p 2 + 

(2.23) 

for an RG shell with and without a Matsubara frequency mode \oj m \ along the n-axis, 
respectively. Once more, Eqs (2.23) and (2.22) coincide when the mode lies exactly along 
the boundary with k — Ak = u> m . 

For (3 — > 0, Eq. (2.8) is again easily recovered. Similarly, for j3 — > oo, uj n can be 
regarded as a continuous variable and we replace the summation by an integral. Dropping 
the subscript (3 in the potential, we obtain 



k 



dk 



S d k d | (3 



2(3 



TV 



t dwJl-^Y 1 lnk + ( fe - Wn ) 2 + ^($) 



+ ln 



S d k d + 
2tt 

S d k d+1 



2% 
S d k d+1 



-/ dx(l-x) lnx 2 + (l-x) + fc ; V o 
cfe (l — x) d 1 + • • • 



+ 



(2.24) 



2dn 



1 + 
1 + 



k 2 

k 2 



'o 
+ • 



In going from the second to the third equation above, we have taken the limit 2ix/(3 — > 0. 
In addition, since the argument in the logarithm does not change appreciably with x over 
the interval, we have dropped its x dependence and moved the expression outside the 
integral. While the expected logarithmic form of the RG flow is regained, the equation 
contains a factor S d /dn instead of the expected S d +i for D = d + 1 from Eq. (2.7). 
Such a constant overall mismatch, however, is harmless and can be eliminated by a simple 
rescaling: k -> c^^'^k and U k ($) -> c 2 /( 2 -^U k ($) with c = S d /S d+1 dn for d ^ 2. 

Up to now we have discussed four types of smearing functions for setting up the finite- 
temperature MRG. In the limit (3 — > 0, all four MRG equations reproduce the expected 
dimensionally reduced result. The agreement can be easily understood from the fact that 
at any arbitrary IR scale k, each of the four different RG shells always contains the static 
n = mode which gives the dominant contribution for small (3, or the high-temperature 
limit. The r ^ modes are strongly suppressed and can be treated perturbatively. In 
this regard, we expect to obtain qualitatively the same physical information from any of 
the four RG equations in the small (3 limit. We shall substantiate this claim later with 
numerical solutions. 
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On the other hand, the asymptotic forms of the four MRG equations in the large (3 
regime all differ from one another. In Figure 6 the comparison between Methods 1 and 3 is 
depicted. The difference may be attributed to their treatment of couplings between |p| and 
the w„^o modes. While the n = sector alone has been shown to be sufficient to attain 
proper scalings for (3 <C 1, to reproduce the less trivial scalings in the (3 ^> 1 regime will 
require a subtle and an intricate account of the feedbacks between |p| and these non-static 
modes. The analyses above show that only Method 3 is constructed in such a way that it 
leads to the desired low-temperature scalings for the running parameters, while rescalings 
or adjustments of the UV cutoff are required for the other methods. 

0.1000 | ■ , ' - tzt^t^- 1 




-7e+04 



Figure 6. Zero-temperature RG flow using Methods 1 and 3. The figure shows that as long as the 
correction to the UV cutoff is taken into account both methods give the same renormalized 
couplings at T=0. 



In passing we comment that there exists a complementary TRG in which (3 acts as 
the running parameter. When (3 is varied infinitesimally, the TRG leads to (see Appendix 
A) 



d(3 



= -/{ilnsi„h(f v /p 2 + C /» p ) 



p 2 + %> + ^v /p2 + c '£p 



coth 



(2.25) 



with k = A — Ak in the |p| integration. This equation is the finite-temperature analog of 
the renown Wegner- Houghton RG [10], and is consistent with the results obtained using 
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the MRG schemes. The general procedure of blocking transformation is not unique, and 
one can certainly construct other types of smearing functions in addition to what have 
been considered. There is an important stipulation, however, for what is considered to be 
an effective RG transformation - that it makes the smallest possible change in going from k 
to k — Ak for any arbitrary k. In other words, the most successful scheme is the one whose 
RG shell contains the least number of modes. Only by doing so is the environment not 
altered grossly. In all four RG prescriptions, the fraction of the modes taken out from each 
blocking step is k = Ak/k. Since each loop is accompanied by a power of k, by taking this 
parameter to be infinitesimally small all higher loop corrections can be suppressed [10], 
[11]. This justifies our MRG equations, which are based on a simple one-loop framework. 
On the other hand, the failure of perturbation theory and the Matsubara shell assumption 
in the high-temperature regime can be readily explained by noting that k ~ 1, which 
implies a dramatic change of the environment. 

One could certainly consider a smearing function such as Pn(uj n ) = Q(n — |n|), where 
the RG evolution is now governed by a difference equation: 



r{™}(^\ , 1 / 1„L,2 , _2 , TT {n}' 



(2.26) 



or 



<>(*) = tfj ra} (*) + ±£ / ln^ + p 2 + ^>"($) 



(2.27) 



In this case eliminating |u> m | implies taking away a fraction of modes, k' = 2/ '(2m + 1), 
from the theory. However, k' is clearly not small at all for small m. Hence this blocking 
method will break down in certain regimes. 



III. NUMERICAL RESULTS 

We now present the numerical results for the four MRG equations discussed in Sec. 
II. By expanding Up t k(&) as a power series in $ 2 

CO 

tW*) = £& m) (*T, (3.i) 

m=l 

followed by a truncation at $ 14 , we transform the non-linear partial differential MRG 
equation into a set of non-linear coupled ordinary differential equations for the couplings, 
{.9 b T /• These equations are then integrated using a fourth-order, adaptive step-size, 
Runge-Kutta integrator. We take d = 3 for definiteness. 

As pointed out in Section II, an entirely continuous RG flow is not obtained when 
|p| and u n are coupled. In particular, the continuous RG flow equations in the form 
of differential equations, e.g. (2.13), can be employed only when the RG shell to be 
eliminated does not span a Matsubara mode. However, when a Matsubara mode, say 
|u; m |, is contained in the shell, a discrete, e.g. (2.15), version must be used. Let the 
distance between k and uj m be 5, where < 5 < Ak. While it is possible to choose 
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the integration bounds in a manner such that the step which spans \oj m \ has 6 = 0, i.e., 
\oj m \ lies exactly on the boundary of the RG shell, it would be interesting to examine the 
dependence of the flow pattern on 6, over which the discrete flow equations are applied. For 
Methods 3 and 4, as 5 — > 0, the discrete equations can be readily shown to reduce to their 
corresponding differential equations. However, Method 2, has a "pathology" due to the 
horizontal sections of the RG shell and is therefore not possible to establish an equivalence 
with its continuous counterpart. In Fig. 7 we give a plot of the temperature dependence 
of the thermal mass parameter as a function of 6 using Method 3. The corrections vanish 
in the limit 5^0. Since having a large 5 will destroy the exactness of the RG flow in the 
same way as having a large k, one must always make 5 as small as possible if not zero. 

5e-05 i ' 1 ' 1 ' 1 ' 1 ' — r^i 




0e+00 1 1 1 1 1 1 1 1 1 1 1 

0.00 0.02 0.04 0.06 0.08 0.10 

T/A 

Figure 7. Dependence of k on 5. using Method 3. 

In Fig. 8 we plot the temperature dependence of the first six coupling constants at 
k = 0. All four schemes show excellent agreement up to oscillations. The oscillations 
obtained using Methods 2, 3, and 4 are due to the boundary conditions on the RG flow. 
For certain T, the constraint u n = A is satisfied exactly at k = A, while for others 

I '"max «/ > 

it is not. Technically speaking one should only compare data points which satisfy this 
constraint because different boundary conditions are mixed in the intermediate points and 
give rise to oscillations in the g^™^ ■ The observed oscillations are a numerical artifact of 
the finiteness of the UV cutoff and signals the elimination of a particular |a> n |. A smooth 
curve may be obtained for methods 3 and 4 if we choose our UV cutoff A to be an integer 
multiple of T and set Ak = 2-7T /M(3 for some integer Af. Irrespective of this technicality, 
it is clear from Fig. 8 that all four methods give qualitatively the same results for the 
temperature dependence of coupling constants. 
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Let us now examine more closely the MRG flow of the relevant coupling constants, 
namely, the mass parameter ^ k and the coupling constant Xp^. At finite temperature, 
according to the prediction of the one-loop perturbative approximation, the running mass 
parameter ftp = fi^ k=0 exhibits a quadratic dependence on T. Applying the four RG 
schemes does not modify its strong T 2 dependence, i.e., incorporating the continuous 
feedback between the modes does not alter the temperature dependence of the thermal 
mass in the high-T limit. 

The situation with the coupling constant A/3 = Xp,k=o, unfortunately, has not been 
so clear. Its high-temperature behavior has been a subject of recent debate [4] [8] . In this 
regime, a naive one-loop calculation breaks down completely and predicts that A/3 decreases 
and turns negative eventually. Including the two-loop contribution does not help either. 
Thus, RG remains the only tool which could shed light on the true behavior of Xp. However, 
the results again differ depending on how RG is implemented and at what scale the coupling 
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constant is denned. In Ref. [8], for instance, it was claimed that when using Method 3, 
A^ decreases at first, but then increases logarithmically. While the author attributed the 
discrepancy between his result and that of Ref. [4] based on Method 1 to the functional 
form of the smearing functions, we find that this is not the case. The actual cause of the 
discrepancy lies in the scale at which the coupling constants are defined, namely, in the 
former, \p^=T was considered, i.e., the coupling constant is defined at the scale k = T, 
whereas Xp t k=o was used in the latter. Had the normalization conditions agreed, the same 
result would have been obtained for all four RG schemes: \(3,k=T decreases, reaches a 
minimum, and then increases logarithmically. However, defining the coupling constant at 
k = shows that \/3,k=o decreases and approximately approaches a constant for large T 
for all four methods. 




Figure 9. Comparison of \p t k with different normalization conditions. 



Certain methods, however, will yield erroneous temperature dependence of the cou- 
pling constants. An example is the Matsubara shell assumption in which a shell of thick- 
ness Ak = 2n//3 is integrated out in each step of the RG evolution. This assumption 
gives k ~ 1 in the high T limit and the method breaks down completely (see Figure 9). 
In addition, this method yields mean field (four-dimensional) critical exponents instead 
of the Wilson-Fisher three-dimensional ones. For example, the exponent f3, defined by 
&p,k ~ \T — T c (k)\P as (3 — > P c {k), where &p,k is the scale-dependent minimum of the 
potential [5], turns out to be (3 = 1.01 ± 0.01 with the Matsubara-shell assumption. 

We comment that the T dependence of the running parameters can also be examined 
from the point of view of the RG trajectories and the fixed-point structure of the theory. 
The presence of the infinite-temperature Gaussian fixed point implies that the RG trajec- 
tory must move in such a way that the fixed-point (/x^*, A£) = (oo, 0) is approached with 
increasing T. However, a careful analytical study reveals that the leading-order increment 
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of jA 2 p is XpT 2 , which means that the flow toward the fixed point must be slowed down as 
Xp continues to decrease according to the coupled relations (see Appendix B): 



X[3 _ Xpup _ Xpfip 
24/3 2 8nf3 2 16tt 2 



Xp = Xr 



3X 2 



sx 2 



16n[ip IQtc 2 



(3.2) 



\47T- 



with ftp = (3 'up. The continuous feedback between ^ and Xp in turn yields a slower 
convergence to their fixed point values. Since for a realistic system the attainment of 
\Tp — > oo is rather unlikely, the corresponding A^ — > limit will not be realized either. 
Therefore, it suffices to approximate Xp by a constant. 

Further argument against the increase of A^ with T can also be made from the consid- 
eration of dimensional reduction, which takes place only if ftp — > in the limit of vanishing 

1/2 

(3 [12]. As we have seen before, the quantity ftp is of order X^ . Hence, if A^ increases 
with T, ftp will also rise and the picture of high-temperature dimensional reduction will 
be destroyed completely. 
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V or C 


littlfc^o Xp cjk T c / fJ,p ci k 


2 


1.054 ± 0.005 


0.527 ± 0.005 


13.79 ± 0.05 


3 


1.171 ± 0.005 


0.585 ± 0.005 


29.35 ± 0.05 


4 


1.345 ± 0.005 


0.672 ± 0.005 


37.23 ± 0.05 


5 


1.504 ± 0.005 


0.744 ± 0.005 


38.78 ± 0.05 


6 


1.352 ± 0.005 


0.675 ± 0.005 


35.50 ± 0.05 


7 


1.312 ± 0.005 


0.655 ± 0.005 


35.20 ± 0.05 



Table 1. Critical exponents and the effective coupling constant as a function of the level of polyno- 
mial truncation. Results for all four smearing functions are the same within the numerical 
accuracy. Errors reported here are due to numerics only and not that resulting from the 
polynomial truncation of the potential. 



Other quantities can be measured to demonstrate the consistency of all four mo- 
mentum blocking schemes. For example, one may examine the case \r R < where the 
symmetry is spontaneously broken. Although Xp vanishes at criticality, the effective mea- 
sure of the interaction is characterized, not by Xp, but by A3 = lim^o liniT^T c Xp^T j \xp^ 
which approaches a constant as T — > T c . The effective coupling strength in the vicinity of 
transition is obtained using the polynomial expansion approximation up to $ 14 for MRG 
and the results are summarized in Table 1. We again find agreement in all four MRG equa- 
tions within their numerical accuracy, and the results are comparable to those reported in 
Ref. [13]. In a similar manner, the critical exponents 7 and v were measured. From Table 
1 we again see that all methods are consistent with one another. However, we caution 
the readers that these results will be less accurate than that obtained without polynomial 
truncation. This is because near criticality field amplitude fluctuations are strong, but 
such truncation applies only to the small fluctuation limit [5] , [14] . 



19 



IV. SUMMARY AND DISCUSSIONS 



In this paper we have constructed finite-temperature RG equations for a scalar field 
theory in S" 1 x R d using a sharp momentum cutoff as the running parameter. We examined 
four different blocking schemes and illustrated how, with an infinitesimally small thickness 
of the RG shell, Ak, the feedback from the higher modes to the lower ones are systemati- 
cally incorporated as they are integrated out. The smallness of the parameter k = Ak/k, 
which represents the fraction of the modes that are eliminated in each blocking step, is the 
most crucial ingredient for the success of RG and constitutes the basis for the "exactness" 
of Eqs. (2.6), (2.14), (2.18) and (2.22), since corrections from higher loops are strongly 
suppressed by additional powers of k. 

The four MRG equations differ essentially in their treatment of the feedbacks between 
modes having the same |p| but different u n . In Method 1 the feedbacks are ignored entirely, 
but the trade-off is that its unconstrained summation over all n leads to a continuous RG. 
On the other hand, Methods 2, 3 and 4 which enforce certain couplings between |p| and 
u n do not lead to a differential function when there are non-vanishing contributions from 
eliminating the modes along the n-axis. However, for Methods 3 and 4, the differential 
forms are recovered by a judicious choice of Ak which puts all these points on the boundary 
of the RG shells. From the numerical results presented in Sec. Ill for the behavior of 
the n-point functions and the critical exponents, one finds remarkable agreement in all 
four methods, thus justifying the use of Method 1. We may therefore conclude that the 
couplings between different frequency modes with the same |p| may be safely ignored, 
provided that k is chosen to be small. Had all the modes been strongly coupled, one 
would have to devise a laborious scheme which eliminates one degree of freedom at a 
time. Nevertheless, it remains an interesting question to see why an "independent-mode" 
approximation for n is sufficient for a small RG shell. 

Knowing that physical properties are independent of the methodology of blocking, 
we find Method 1 to be most advantageous for numerical computation since it yields a 
continuous flow, with no oscillatory behavior in the high T limit. On the other hand, to 
examine the scaling behavior of the running parameters in the low T regime, Method 3 
would be more convenient since it naturally leads back to the zero-temperature (d + 1)- 
dimensional theory. 

Using the MRG (and TRG), we have examined the high-temperature behavior of 
coupling constant Xp. Based on the fixed-point structure of the theory and the numerical 
results, it is clear that \p,h=o must decrease with T but can be approximated by a positive 
constant in the limit flp <C 1. The discrepancy in the literature concerning the behavior of 
Xp ; k can be reconciled with a precise renormalization condition (viz, k = or k = T). We 
emphasize again that integrating out entire Matsubara shells at each RG step is ill-advised 
since it yields k ~ 1 in the large T limit and the exactness of the RG flow is lost. 

Various interesting directions can be pursued with our RG equations: Numerous crit- 
ical exponents have been measured by Method 1 with the inclusion of wavefunction renor- 
malization Z^ fc($) in [5]. However, the accuracy in the polynomial expansion can certainly 
be improved. Derivations of the RG equation for Zp^{®) using Methods 2, 3 and 4 can 
also provide a consistency check with Method 1 [15]. It would also be interesting to 
compare our predictions for the universal amplitude ratios with those obtained using the e 
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expansion and the lattice methods. The one-component RG equation Eq. (2.6) can readily 
be generalized to N components: 

k dU^ = _^{ (iV _ l) ^ k , + u , k w + 2 ln [! _ W^)]) 

+ /3 v /p + ^ fc ($) + ^ fc ($)$2 + 21n 

(4.1) 

where the prime notation now denotes the differentiation with respect to $ 2 /2 = $ a $ a /2 
with a = l,---,iV. Critical exponents corresponding to the XY (N = 2) model and 
Heisenberg (N = 3) model can now be extracted in a similar manner. In addition, the 
equation may also provide information on the interplay between the longitudinal and the 
transverse Goldstone modes. One may also explore different IR scaling regimes and search 
for the possible non-universal behavior which may arise due to the emergence of new 
relevant operators that are suppressed in the UV scaling regions [16]. Nonperturbative 
modification of the scaling laws due to the presence of a condesate can also be addressed. 
Gauge theories, too, can also be further investigated [17] [18]. These issues will be 
presented in forthcoming publications. 

ACKNOWLEDGMENTS 

S.-B. L. is grateful to J. Polonyi, V. Branchina, J. Alexandre and A. Bodor for stimu- 
lating discussions. M.S. would like to thank D. Furnstahl and E. Braaten for discussions as 
well. This work is supported in part by funds provided by the National Science Council of 
Taiwan under contract #NSC-87-2112-M-194-004, and by the National Science Foundation 
under Grants No. PHY-9511923 and PHY-9258270. 

APPENDIX A. TEMPERATURE RENORMALIZATION GROUP 

In Sec II, the MRG prescription was explored from the point of view of momentum 
space construction using an IR cutoff k as the running parameter. The external variable, 
temperature, has been taken as an input parameter. To examine the behavior of a physi- 
cal system at finite temperature, there exists an alternative temperature renormalization 
group (TRG), in which f3, the inverse of temperature, plays the role of the flow param- 
eter [19]. TRG provides a direct investigation of the response of a system when T is 
varied. Quantities at different T and fixed k can also be related to each other by this finite 
renormalization procedure. 

To construct the TRG flow equation, we again begin with the perturbative one-loop 
expression for the blocked potential given in Eq. (2.3). Differentiating the expression with 
respect to (3 yields 
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By varying (3 infinitesimally from /?—>/? — A[3, i.e., the system is heated up gradually 
until the desired temperature scale is reached, we arrive at the following improved TRG 
equation: 







dU, 
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(3 ln Sm H 2 , 



P 2 + Ufa+(3 



d_ 

op 



p 2 + "I* 



coth 



(A.2) 



By assuming a small RG shell for the |p| integration, or equivalently, k = A — A/c, higher 
loop corrections again are suppressed. 

One can readily verify the equivalence between the TRG and the MRG formalisms by 
observing that the same expression 
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(A.3) 



is obtained by differentiating (A.2) with respect to ln/c or (2.6) with respect to In/?. How- 
ever, since Eq. (A.2) is an integro-differential equation, direct numerical solution is more 
cumbersome than the MRG equations. Notice that for large (3, the results derived from 
(A.2) and (A.l) should not differ appreciably. However, in the small (3 regime, Eq. (A.l), 
which is derived using the independent-mode approximation, fails completely. 

The (3 dependence of the running parameters in the coupling constant space can be 
examined by differentiating (A.2) with respect to $ and setting $ = 0. To 0(Xp ! k), we 
have 



p 2 + ^ 
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(A.4) 
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(A.5) 

where C.T. denotes the appropriate counterterms needed to subtract off the divergences 
in the momentum integration. Note that these counterterms are necessarily temperature- 
dependent [20] . 

Unfortunately, the coupled integro-differential equations are still too difficult to be 
solved exactly; even analytical approximations are restricted only to certain regimes. How- 
ever, in the high-temperature regime where (3 is small, we ignore the p dependence in the 
mass and coupling constants and obtain 
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and 
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for vanishing k. The above two RG equations will allow us to examine the fixed-point 
structure of the theory, as we shall see in Appendix B. 



23 



APPENDIX B. RENORMALIZATION GROUP AND THE FIXED POINTS 



The complicated RG flow can be greatly simplified in the neighborhood of fixed points 
where a linearized RG prescription becomes possible and only the relevant operators need 
to be retained. The irrelevant operators can be safely ignored. The classification of rel- 
evancy, however, depends on the fixed point; what is relevant around a particular fixed 
point may turn irrelevant as it flows into another one. Thus, the linearized prescription 
may be used only when the system is sufficiently close to a fixed point. To describe the flow 
further away from the fixed point, a global nonlinear RG analysis is generally required. 
Once all the fixed points are located, the general RG flow of the theory can be mapped 
out. Below we examine the fixed-point structure of both MRG and TRG. 

1. MRG Approach 

For MRG, we shall work with Method 1 since its RG evolution is continuous. It is 
convenient to employ the dimensionless quantities: 

Up, k @) = (3k~ d Up^), $ = /? 1 /2/ c -(d-2)/2 $) 

^(*) = = ^- m/2 k- d+ ^ d -^U^^), (B.l) 

fiik = U™(0) = k- 2 ^, Xp, k = C#j>(0) = p-'k^X^ , 



where, for the purpose of the present paper, we set the anomalous dimension rj to be zero. 
Eq. (2.6) can now be written as 



kd k - -(d-2)<f>d$ + d 




'l + U% k + 21n[l - e - V 1+0 '^] | , (B.2) 

which implies the following fixed-point equation: 

-\{d- 2)*U}[ k + dUl k = + % + 21n[l - e^V^Z] j , (B.3) 

with Up k (&) being the fixed-point potential. However, by making a polynomial expansion 
of the blocked potential Up ik (^) and keeping track of only the evolution of the leading- 
order Jl 2 p k and Xp >k , it is possible to identify the fixed points on the (ftp k , Xp jk ) plane. 
Consistency with reflection symmetry around the origin demands that only even-power 
terms are present and we have 

oo -(2m) 

°f>>*W = E §^ 2 ™ 4 2 r } = ^ 2 T\0) , (B.4) 

m=l ^ ' 

with 9pl = P>}, k and and gf^ = Xp, k . 
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In the high-temperature limit where C 1, the manifold changes from 5" 1 x R d to 
R d : exhibiting a critical behavior governed by the <i-dimensional fixed point. Thus, one 
expects 

( S 1 xR d — ► R d 
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Concentrating only on the running of /i| fc and the evolution equations can be ap- 

proximated as 
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Eq. (B.7) then simplifies to 



fc* -(l-f)(2f+^) 



k dk~ 



Ar 



5S d 



)■ 



(B.9) 



where e = 4 — d. Fixed points can now be located. There are two trivial fixed points 
of Gaussian nature located at (f*,£*) = (0,0) and (1,0). The former is the usual trivial 
"finite" Gaussian fixed point at which all running parameters vanish. On the other hand, 
the latter describes an infinite (temperature) Gaussian fixed point which can only be 
reached when T — > 00. In addition, for D = d < 4 or e > 0, we have the IR Wilson-Fisher 
(WF) fixed point (f*,£*) = (—e/6,2e/3S d ) [21], which coincides with the usual trivial 
Gaussian fixed point for e = 0. Thus, we have £* = 16iv 2 e/3 for d = 4. While Eq. (B.9) 
is already diagonal around the infinite Gaussian fixed point, it is possible to introduce the 
new variables [22] 

x = ^£, y = f+ n J d J , (B.10) 



and transform (B.9) into 
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The new expression is now diagonalized around the finite Gaussian fixed point (0, 0), with 
the infinite Gaussian and the WF fixed points located, respectively, at (0, 1) and (1, 0) on 
the (x, y) plane. 

On the other hand, in the low T limit, one expects 
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Here we see that the only fixed point is the trivial finite Gaussian (/i|* d+1 , A£ d+1 ) = (0, 0), 
for d > 3. While the finite Gaussian fixed point can be accessed in both low and high T 
regimes, the WF and the infinite Gaussian fixed points are reachable only in the large T 
limit. The rationale for the inability to detect these two fixed points at low T is that they 
are characteristic of R d with d < 4, and yet the effective manifold is R d+1 for small T. 

2. TRG Approach 



The fixed-point structure associated with the flow of TRG can be examined in a similar 
manner. We focus again only on the evolution of the leading-order running coupling 
constants fip and A^ at k = 0. To simplify the task of solving the coupled integro- 
differential equations Eqs. (A. 6) and (A. 7), we choose the manifold to be S 1 x R s for 
definiteness. 

In the limit (3fip <C 1, the RG equations can be approximated as 
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and [23] 
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yV + 1 e ^v / y 2 + 1 - l 



3A| 



2tt 4tt 2 
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Note that /3 is now the running parameter of the RG. One may actually derive the above 
evolution equations by differentiating the coupled self-consistent equations 



2 2,^/3 



24/3 2 8tt/3 2 16tt 2 



\r = A 



3A 2 



3A; 



1677//^ 167T 2 



In 



+ 7 
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which are valid for a small but positive \ip [5]. 

Upon substituting the familiar leading-order result 







A/3 



+ 



or 



2 2 



into Eq. (B.16), the fixed points may now be located by solving 

>9Xp 



(B.20) 



= fi- 



d/3 



3A 2 , A 
16tt 2 



( X 24^ 2 )( 1 Up) 



(B.21) 



Since the equation is derived for p,p <C 1, we require ftp < tt. Two fixed points may now 
be located. The first one is the usual trivial finite Gaussian fixed point (fip, \p) = (0,0). 
As for the second one, we see that for A/3/24/2 2 — > 1, or equivalently f3 — > 0, the theory 
approaches the infinite Gaussian fixed point (fif3,\p) = (oo,0). Its Gaussian nature is 
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1/2 

readily seen by noting that up ~ A^ , and thus, when \xp — > oo and p — >■ such that 
/x^ — > 0, A^ must be vanishingly small. 

In general one may regard phase transitions as a "high-temperature" phenomenon 
since \ip — > in the vicinity of criticality. Although we are able to detect two Gaussian 
fixed points, the important WF fixed point which characterizes the critical behavior of 
the system seems to be missing. A crucial observation here is that in the large T regime 
one expects the effective degrees of freedom to crossover from four-dimensional to three- 
dimensional, and hence Xp is not an accurate measure of the effective interaction strength 
of the theory. Instead, the interaction is accounted for by Xp^ = XpT, or Xp^ = Xp/ftp in 
the dimensionless form. Thus, in terms of ftp and Xp,3, we have 



R dfl/3 _ X/3,3 

,d\fji T / A 



(B.22) 



The much richer structure associated with the flow of A^ now reveals the existence of 
the additional nontrivial fixed point. Besides the usual finite (0, 0) and infinite Gaussian 
(1,0), we have, keeping only the leading-order contribution in the bracket, 

(^ 3 ) = (f , ^ ), (B-23) 

The numerical value for A^ 3 can be compared with that tabulated in [13]. Therefore, 
contrary to the well known notion of triviality for T = at dimensionality D = 4, the 
theory in the high T limit is seen to crossover to an interacting theory characterized by a 
nontrivial fixed point of D = 3. 

We comment that in the TRG formalism, by first taking the limit k — > in Eqs. 
(A. 4) and (A. 5), we are left with one arbitrary parameter /?, the circumference of the 
submanifold S 1 , with all other dimensions being infinite. Since any TRG transformation 
now formally corresponds to changing the value of /?, it is natural to expect the extrema 
(3 = and P = oo to be fixed points that characterize, respectively, the three-dimensional 
and four-dimensional limits of the theory. While D = 4 for p,p ^> 1, it becomes three 
in the limit jlp 1. Therefore, the dimensional crossover scale in this TRG formalism 
can be given qualitatively as p,p ~ 1. This implies that when ftp falls below unity new 
effective three-dimensional degrees of freedom need to be employed in order to give a 
correct account of the theory. However, had the new degrees of freedom associated with 
D = 3 not been used, the WF fixed point would be missed. In fact, the success of a given 
RG scheme in characterizing the crossover phenomenon is attributed to its capability of 
tracking the effective degrees of freedom. We also note that the method of resumming 
daisy and super daisy diagrams proposed by Dolan and Jackiw [24] is similar to TRG; 
however, the important distinction is that in the former approach, the coupling constant 
has always been fixed at its zero-temperature renormalized value and consequently, the 
prescription only gives the two Gaussian fixed points and yields the mean-field exponents. 
Thus, we see that only when the running of the effective coupling constant \p^ is taken 
into account can one access the nontrivial WF fixed point. 
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FIGURE CAPTIONS 



Figure 1 
Figure 2 

Figure 3 
Figure 4 
Figure 5 
Figure 6 

Figure 7 
Figure 8 

Figure 9 
Table 1 



Schematic diagram of blocking Method 1. 

Schematic diagram of blocking Method 2. The thickness of the RG shell in the 
n-axis is zero for Ak — > (exaggerated here), but a Matsubara mode is eliminated 
after M continuous iterations where [A/"/3A/c/27r] = 1. 

Schematic diagram of blocking Method 3. 

Plot of g^ifi) for Method 3. 

Schematic diagram of blocking Method 4. 

Zero-temperature RG flow using Methods 1 and 3. The figure shows that as long 
as the correction to the UV cutoff is taken into account both methods give the 
same renormalized couplings at T = 0. 

Dependence of ^ h on 8. using Method 3. 

Temperature dependence of symmetric phase for m = 1, •••6. The 

solid, dotted, dashed lines, and dot-dashed lines are for Methods 1, 2, 3, and 4 
respectively. 

Comparison of Xp^ with different normalization conditions. 



TABLE 

Critical exponents and the effective coupling constant as a function of the level 
of polynomial truncation. Results for all four smearing functions are the same 
within the numerical accuracy. Errors reported here are due to numerics only 
and not that resulting from the polynomial truncation of the potential. 
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